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Note on a standard program for calculation of one-electron 
molecular integrals of overlap type by use of the Swedish 
electronic computor BESK’ 


By Sric¢ FLopMARK 


SUMMARY 


All the combinations of overlap integrals between the hydrogen-like wave functions of ls, 
.2s, 2p, 3s, 3p and 3d types can be calculated by the use of this program. This is specially constructed 
for overlap integrals between equal atoms, but can be used even for unequal atoms, except when 
the orbitals of the different atoms have equal quantum numbers. The atomic radius and the effec- 
tive charges of all the orbitals must be put on the parameter tape together with the number of 
those integrals which will be calculated. The computation time is about 3 secs per integral if 
|B| > 2 or if B = 0, where f is given in (8). If || < 2 the time will be about 10-15 secs per integral. 
A special arrangement has been made in order to vary the inter-atomic distance. Partial results 
as the values of An (a) and Bn (f) for n = 0,1,...,8, which are defined in (3), can also be given 
by the program, if a special code notation will be put on the parameter tape just before the number 
of the integral. 


The method of calculation 


The wave functions ‘’ are numbered by an index K = 0, 1, ..., 9 as follows: 
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1 
W ls 2s 2px 2p0 33 3px 3po 3dd 3dx 3do/ a 
The overlap integrals are numbered by two indices KL. The z- and o-directions of 
the two participating atoms are chosen to be equal (parallel coordinate systems). 
The direction of o goes from the first to the second atom. 
All the integrals KZ can be written as a linear combination of expressions of the 
following type: 
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The integers C,, v, and y, of a single integral are calculated by BESK from the 

98 at less than 10-* secs. All the ¢7,"’s are placed on the program tape and will be 
cht over to the magnetic drum of BESK before the calculations start. The actual 
set of ¢ has been chosen by BESK according to the number KL of the integral on 
the parameter tape. 


1 A new program giving all types of molecular one-electron two-center integrals of the K, L 
and M shells will be finished in a short time. 
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The calculations of the functions 


A, (a) = 


om, 8 


pertas,| 


: (3) 
By (B)= J ne P" dn | 
“4 
can be made according to the well-known formulae («+0), 
A, (0) =~ [vy Aya (a)+e™*] Ay (a) = > (4) 
a om 

By, (B) = —Au(8)—(—1) 4u(—B) (89), (5). 

as . 

By (0) = (6) 


These functions will be computed by BESK for yw, » = 0, 1, ... 8 within a fraction 
of a second. However, the possible accuracy of the exponential functions at present 
at BESK, which corresponds to 12 decimals, will not at all be sufficient for calcula- 
tion of the B-functions according to (5), especially for large u-values and small 
B-values. It might happen that even the order of magnitude will be incorrect. To- 
correct this, it is possible to use the following expansion of Bin series, where cancelling 
terms in the ordinary Taylor series are omitted: 


vf 
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Working with (7) it might, however, take 30 secs or more for BESK to calculate 
an integral value. Another inconvenience is the slow convergence for large values 
of 8, which means many additions and low accuracy. In the present program (7) 
will be used when |f|<2 and (5) when |8|> 2. This combination always gives 8. 
correct decimals. 

Before the calculations of the separate integrals, the values of the parameters «. 
and # and the normalizing factors N will be calculated for all possible combinations. 


KL according to 
Zk Z 
xe = R (2 + =) ? 


NK NL 
Ber =R (==- “| (8). 
Ne NL 


Nei =Crz(RZx)"k** (RZ_) "2*4, 


where half the inter-atomic distance R and the effective charges Z must have been 
put on the parameter tape. The principal quantum numbers n and the constants 0” 
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are chosen according to the KL-number on the parameter tape. This number will 
also make it possible for BESK to choose the actual linear combination of (2). 

The resulting value will be type-written or punched on a tape, which takes about 
one second. If the new rapid punch is in order, this time will be almost entirely saved. 
The total time per integral, printing included, will thus be about 3 sees if |f| > 2 
(or if 6 = 0) and up to about 15 secs if |8|<2. The rapid punch saves 1 sec per inte- 
gral. 

The program will be kept available at BESK and can be used by everybody. The 
price depends on the time used on BESK and will at present be 6 Sw. kronor a 
minute. Except for the calculating time, there will be about 2 minutes to get the 
program into the machine. Details concerning the managing of the program can be 
sent to those interested. 
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